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ABSTRACT 



■ We report NIR broad and narrow band photometric observations in the direction 

■ of the IRAS17149-3916 source that reveal the presence of a young cluster of massive 
stars embedded in an HII region coincident with RCW121. These observations, together 
with published radio data, MSX and Spitzer images were used to determine some of the 



o 

(N 



, physical parameters of the region. We found 96 cluster member candidates in an area 

of about 1.5 X 2.0 square arcmin, 30% of them showing excess emission in the NIR. 
IRS 1, the strongest source in the cluster with an estimated spectral type of 05V-06V 
ZAMS based on the color-magnitude diagram, is probably the main ionizing source of 
' the HII region detected at radio wavelengths. Using the integrated Br7 and the 5 GHz 

OO . flux densities, we derived a mean visual extinction Ay = 5.49ibf'32 magnitudes. From 

. the observed size of the Br7 extended emission, we calculated the emission measure 

, E = 4.5 X lO'^^cm"^ and the electron density He = 2.6 x lO^cm"^, characteristic of 

Q I compact HII regions. 

Qh' Subject headings: stars : formation - stars: pre-main sequence - infrared : stars - ISM: 

O ' HII regions - ISM: dust, extintion 

w 

1. Introduction 

X 

' Young massive stars are found inside dense molecular clouds generally suffering high extinction 

at optical wavelengths, mainly due to the large amount of gas and dust remaining from the star 
forming process (Churchwell 1990; Cesaroni 1991). Color-color and color-magnitude diagrams 
from near infrared photometry (NIR) allow the determination of the reddening of each individual 
star (at least for those that do not present intrinsec excess emission in the NIR). Longslit i^-band 
spectra have shown that a large number of those NIR sources present broad Br7 emission and no 
photospheric features (Armand et al. 1996; Hanson, Luhman & Rieke 2002), typical of very young 
pre-main sequence stars surrounded by disks and/or cocoons (Churchwell 1990). 

On the other hand, narrow band imaging observations of star forming regions (e.g. centered 
in the Br7 line and in the adjacent continuum) are powerful tools to study the emission from the 
extended surroundings. In fact, by comparing the Br7 flux density with the free-free radio emission. 
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it is possible to estimate the absorption (Hanson, Luhman & Rieke 2002). Moreover, using the 
measured size of the Br7 emitting region and the low resolution radio flux density, it is possible to 
obtain the emission measure and mean value for the electron density (Mezger & Henderson 1967). 
We have to take into account that the measured size of the Br7 emitting region is limited by the 
sensitivity of the NIR detector systems. The recent near to mid infrared images obtained by the 
Spitzer space telescope ^, represent an additional constrain to the size of the ionized region, since 
they delineate the dusty border of the HH region. 

RCW121 located at 1=348?2 and b=-l?0 (Rodgers, Campbell & Whiteoak 1960) is a strong 
(11.9 Jy) and compact (about 1 arcmin angular diameter) 5 GHz continuum source (PMNJ 1718- 
3916), observed in the H109a and HllOa recombination lines by Caswell & Haynes (1987). The 
CS(2-1) transition, typical of dense molecular clouds, was also detected by Bronfman, Nyman &: 
May (1996) and methanol maser emission at 6.7 GHz, generally associated to high mass stellar 
forming regions, was observed with the Parkes radio telescope by Walsh et al. (1997). Maser 
emission from excited OH in the ^Hi/2 J=l/2 state was found by Cohen, Masheder k. Caswell 
(1995) at all three transitions (4765, 4750 and 4660 MHz). Dutra ct al. (2003), by visual inspection 
of 2MASS images, found a cluster of stars in the direction of RCW121; this fact, together with the 
presence of an IRAS source with colors of compact HII regions (Wood & Churchwell 1989) indicates 
that star formation has recently occurred. Despite of the overwhelming evidences pointing to the 
existence of young and massive stellar objects associated with RCW121, no detailed studies of the 
stellar content of the region were made up to the present. 

This work is part of a larger survey (Roman-Lopes et al. 2003; Roman-Lopes & Abraham 
2004a,b) , aimed at the detection of the exciting stars of compact HH regions associated with IRAS 
point sources which are also strong CS emiters (Bronfman, Nyman &; May 1996). In the next 

sections we present the results of NIR broad and narrow band photometric observations of a young 
cluster of massive stars found in the direction of the IRAS17149-3916 source. In Section 2 we 
describe the NIR observations; in Section 3 we present the results from the broad and narrow band 
NIR filters and in Section 4 we summarize the work. 

2. Observations and data reduction 

The NIR imaging photometric observations were performed in 2003 May in the direction of 
IRAS point source 17149-3916, using the near infrared camera of the Laboratorio National de 
Astrofisica, Brazil, mounted on the 0.6 m Boiler &; Chivcns (BC) telescope: the resulting plate 
scale of the images is 0.48"pixeri. The images were made using the J and H broad band filters 
and the K narrow band filter Ci, centered at 2.14^m continuum. Since meaningful information 
about the ionized regions can be obtained from narrow band Br7 images (Roman-Lopes & Abraham 
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2005), we also made observation in 2005 May, using the near infrared camera, mounted on the 1.6 
m Perkin-Elmer (PE) telescope; the resulting plate scale of the images is 0.24"pixel~^. 

The observations (mean FWHM values of the point-spread functions, the total integration 
times, etc ) arc summarized in Table 1. The NIR images were reduced using programs within the 
IRAF image reduction package system. Each image was corrected for nonlinear response, dark- 
subtracted, flattened, and masked to eliminate bad pixels. We observed calibration stars from the 
list of 2MASS standards at similar air mass of our science targets. In addition we checked our 
photometric results by comparing it with the 2MASS photometry for common bright field stars. 
Details about the calibration and reduction procedures can be found in Roman-Lopes et al. (2003). 

3. Results 

Figure 1 shows the J, H and K near infrared images obtained during the 2003 observing run 
(upper, central and lower panels, respectively). 

3.1. The stellar cluster 

To separate the cluster population from the field stars we analyzed the stellar colors in two 
distinct areas as can be seen in Figure 1: one, labeled "cluster" and measuring about 2.2' x 2.4', 
encloses the IR nebula; the other, labeled "control" and measuring about 8.8' x 1.3', includes the 
eastern and southern parts of the images. From the [J — H) versus {H — K) comparative diagrams, 
shown in Figure 2, we can see that a significant number of "cluster" sources have excess emission 
at the longer wavelengths, probably due to the presence of warm circumstellar dust or disks (Lada 
& Adams 1992). On the other hand, the majority of the sources in the "control" diagram lie inside 
the reddening band for the Galactic field stars. They can be divided in two groups separated by a 
gap around (J — H) ~ 1.4: the first is formed by foreground field stars while the other is composed 
probably by giants stars located behind the molecular cloud. 

Table 1. Summary of the NIR observations 



Filter Date Telescope ITIME FWHM(") 



J 


May 2003 


BC 2 


840 


1.9 


H 


May 2003 


BC 


420 


2.0 


K(C1) 


May 2003 


BC 


7200 


2.1 




May 2005 


PE 3 


360 


1.7 


K{Cl) 


May 2005 


PE 


360 


1.7 



^ Boiler & Chivens 0.6 m telescope 
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We assumed as candidate cluster members all NIR sources that lie along the early type redden- 
ing band and those showing excess emission in the "cluster" color-color diagram. For the infrared 
sources not detected in the J band a different criteria had to be used. In order to distinguish the 
genuine candidate cluster members from the field stars, we analyzed for the two selected areas, the 
number of stars N{H — K) with {H — K) colors as a function of {H — K). The results with the 
counts normalized to the cluster region area are shown in Figure 3; we can see in the cluster area 
an excess of counts in the color ranges {H — K) ~ — 0.75 and {H — K) > 1 . The first range 
corresponds to the cluster sources and the second group can be formed either by genuine cluster 
members or by giant stars located behind the molecular cloud. For this reason, only six stars not 
detected in the J band, with colors {H — K) ~ — 0.75 or {H — K) > 1 and located inside the 
nebula region were taken as candidates cluster members. 

In Figure 4, we show the K band contour diagram in which the relative position of each 
selected source is indicated. There, we show the objects that present excess emission in the NIR 
(filled triangles), those without (open triangles), as well as the sources detected only in the H and 
K bands (filled squares). 

We assumed a kinematic distance to the cluster of 1.6 kpc (Walsh et al. 1997), for a galacto- 
centric distance of 8.5 kpc. Using this value, the luminosity of the stars in Table 2 can be inferred 
from the J versus (J — H) color-magnitude diagram shown in Figure 5. The position of the main- 
sequence is plotted together with each spectral type earlier than F8V. The intrinsic colors were 
taken from Koornneef (1983) and the absolute J magnitudes were computed from the absolute 
visual luminosity for ZAMS given by Hanson et al. (1997). The reddening vector for a ZAMS BO 
V star, as in Fitzpatrick (1999), is shown by the dashed line together with the positions of visual 
extinction Ay = 10, 20 and 30 magnitudes indicated by filled circles. 

From the color-magnitude diagram, we can see that IRSl, IRS5, IRS9 and IRS12 are the 
most luminous objects without excess emission in the whole region, being the best ionizing source 
candidates. An inspection of the relative position of these sources shows that only IRSl is located 
in the central part of the infrared nebula; the others lie at the border of the extended emission 
region. In this sense IRSl, an 05V-06V ZAMS from Figure 5, probably is the main source of 
Lyman continuum photons that are ionizing the HII region detected at radio wavelengths. 
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Tablc 2. List of the selected near-infrared soTirces 



IRS a(J2000) (5(J2000) JcamlV HcamlV KcamlV 



1 


17: 


18: 


:25, 


.51 


-39: 


:19; 


;08. 


.8 


8, 


.93 


±0. 


,02 


8. 


.36 


±0. 


.02 


8, 


.02 


±0. 


.02 


2 


17: 


18: 


;25, 


.73 


-39: 


:18: 


;27, 


.0 


11, 


.81 


±0. 


,02 


9, 


.64 


±0, 


.02 


8 


.45 


±0, 


.02 


3 


17: 


18: 


:22, 


.27 


-39: 


:18: 


:42, 


,5 


11, 


.60 


±0. 


,02 


10, 


,23 


±0, 


.02 


9 


.06 


±0, 


,02 


4 


17: 


18: 


:22, 


.92 


-39: 


:18: 


:22, 


,7 


15 


.03 


±0, 


,03 


11, 


,70 


±0, 


,02 


9 


.72 


±0, 


,02 


5 


17: 


18: 


:20, 


.59 


-39: 


:18: 


:07. 


.2 


11, 


.87 


±0. 


,03 


10, 


.66 


±0. 


.02 


10, 


.06 


±0. 


.02 


6 


17: 


18: 


:30, 


.72 


-39: 


:18; 


;40. 


.0 


11, 


.43 


±0. 


,02 


10. 


.75 


±0. 


.02 


10 


.26 


±0. 


.02 


7 


17: 


18: 


:25, 


.66 


-39: 


:19; 


;09. 


.8 


11, 


.61 


±0. 


,03 


10. 


.86 


±0. 


.03 


10 


.32 


±0. 


.04 


8 


17: 


18: 


:26 


.62 


-39: 


:17; 


;51, 


.7 


12 


.05 


±0. 


,02 


10. 


.98 


±0, 


.03 


10 


.40 


±0, 


.02 


9 


17: 


18: 


:25, 


.34 


-39: 


:18: 


;24, 


,8 


13 


.53 


±0. 


,03 


11, 


,47 


±0, 


,02 


10 


.42 


±0, 


,02 


10 


17: 


18: 


:25, 


.32 


-39: 


19: 


:09, 


.5 


11, 


.87 


±0. 


,03 


11, 


,35 


±0, 


,02 


10 


.60 


±0, 


,08 


11 


17: 


18: 


:28, 


.75 


-39: 


19: 


:15, 


,6 


11, 


.23 


±0. 


,02 


10, 


,96 


±0, 


,02 


10 


.69 


±0, 


,03 


12 


17: 


18: 


:30, 


.36 


-39: 


19: 


:53. 


.0 


14, 


.97 


±0. 


,04 


12. 


.16 


±0. 


.02 


10, 


.75 


±0. 


.03 


13 


17: 


18: 


:30, 


.70 


-39: 


:19; 


;52. 


.7 


11, 


.52 


±0. 


,03 


11. 


.12 


±0. 


.02 


10 


.86 


±0. 


.03 


14 


17: 


18: 


:25, 


.39 


-39: 


:19; 


;12. 


.0 


12, 


.30 


±0. 


,03 


11. 


.36 


±0. 


.03 


10 


.91 


±0. 


.03 


15 


17: 


18: 


:31 


.90 


-39: 


19; 


;40, 


.4 


12 


.27 


±0. 


,03 


11, 


.33 


±0, 


.03 


10, 


.93 


±0. 


.03 


16 


17: 


18: 


:32 


.04 


-39: 


:17: 


:54, 


,2 


13, 


.87 


±0. 


,03 


12, 


,06 


±0, 


,03 


11 


.08 


±0, 


.03 


17 


17: 


18: 


:31, 


.68 


-39: 


18: 


:24, 


,1 


12, 


.28 


±0. 


,03 


11, 


,48 


±0, 


,03 


11 


.13 


±0, 


,03 


18 


17: 


18: 


:25, 


.18 


-39: 


:18: 


:46, 


,8 


14, 


.87 


±0. 


,03 


12, 


,44 


±0, 


,03 


11 


.20 


±0, 


,03 


19 


17: 


18: 


;24, 


.48 


-39: 


:19: 


;39. 


.7 


11, 


.70 


±0. 


,03 


11. 


.43 


±0. 


.03 


11, 


.22 


±0. 


.03 


20 


17: 


18: 


:22, 


.42 


-39: 


:18; 


;21. 


.2 


12, 


.77 


±0. 


,03 


11. 


.74 


±0. 


.03 


11 


.21 


±0. 


.03 


21 


17: 


18: 


:32, 


.98 


-39: 


:19; 


;32. 


.2 


12, 


.04 


±0. 


,03 


11. 


.58 


±0. 


.03 


11 


.27 


±0. 


.03 


22 


17: 


18: 


:25, 


.75 


-39: 


19; 


;55, 


.9 


14, 


.35 


±0. 


,03 


12, 


.33 


±0, 


.03 


11 


.41 


±0, 


.03 


23 


17: 


18: 


:21, 


.70 


-39: 


19: 


:15, 


,2 


14, 


.78 


±0, 


,03 


12, 


,66 


±0, 


,03 


11 


.45 


±0, 


,03 


24 


17: 


18: 


:21, 


.60 


-39: 


:18; 


:31, 


.7 


14, 


.91 


±0. 


,03 


12. 


.63 


±0, 


.03 


11 


.46 


±0, 


.03 


25 


17: 


18: 


:29, 


.40 


-39: 


:17; 


;56. 


.0 


12, 


.30 


±0. 


.03 


11. 


.78 


±0. 


.03 


11 


.47 


±0, 


.03 


26 


17: 


18: 


:30 


.91 


-39: 


:18; 


;33. 


.8 


12, 


.66 


±0. 


.03 


11. 


.89 


±0. 


.03 


11 


.49 


±0, 


.03 


27 


17: 


18: 


:20 


.59 


-39: 


:18; 


;13. 


.0 


12 


.02 


±0. 


,03 


11. 


.72 


±0, 


.03 


11 


.50 


±0, 


.03 


28 


17: 


18: 


:22, 


.44 


-39: 


19: 


:20, 


,6 


12, 


.89 


±0, 


,03 


12, 


,07 


±0, 


.03 


11 


.52 


±0, 


.03 


29 


17: 


18: 


:25, 


.97 


-39: 


19: 


:26, 


,8 


13, 


.48 


±0, 


,03 


12, 


,29 


±0, 


.03 


11 


.59 


±0, 


.03 


30 


17: 


18: 


:26, 


.16 


-39: 


:17: 


:42, 


,7 


15, 


.25 


±0, 


,03 


12, 


,88 


±0, 


.03 


11 


.61 


±0, 


.04 


31 


17: 


18: 


:23, 


.33 


-39: 


19; 


:08. 


.0 


13, 


.45 


±0. 


.03 


12. 


.35 


±0, 


.03 


11, 


.76 


±0, 


.03 


32 


17: 


18: 


:31, 


.37 


-39: 


:18; 


;10. 


.4 


14, 


.98 


±0. 


.03 


12. 


.84 


±0. 


.03 


11 


.84 


±0. 


.03 


33 


17: 


18: 


:23, 


.83 


-39: 


:19; 


;04. 


.4 


12, 


.98 


±0. 


.03 


12. 


.31 


±0, 


.03 


11 


.86 


±0. 


.03 


34 


17: 


18: 


:29, 


.69 


-39: 


19; 


;02, 


.3 


12 


.74 


±0. 


,03 


12, 


.23 


±0, 


.03 


11 


.90 


±0, 


.03 


35 


17: 


18: 


:25, 


.54 


-39: 


19: 


:23, 


,2 


13, 


.55 


±0, 


,03 


12, 


,48 


±0, 


.03 


11 


.90 


±0, 


.03 


36 


17: 


18: 


:29, 


.81 


-39: 


19: 


:07, 


,3 


15, 


.62 


±0, 


,04 


13, 


,21 


±0, 


.03 


11 


.91 


±0, 


.03 


37 


17: 


18: 


:21, 


.48 


-39: 


19: 


:13, 


,8 


13, 


.05 


±0, 


,03 


12, 


,39 


±0, 


.03 


11 


.92 


±0, 


.03 


38 


17: 


18: 


;27, 


.89 


-39: 


:17: 


:58. 


.6 


15, 


.28 


±0. 


.04 


13. 


.02 


±0, 


.03 


11, 


.93 


±0, 


.03 


39 


17: 


18: 


:21, 


.77 


-39: 


:19; 


;22. 


.1 


15, 


.36 


±0. 


.05 


13. 


.24 


±0. 


.03 


11 


.93 


±0, 


.03 


40 


17: 


18: 


:24, 


.48 


-39: 


:19; 


;21. 


.4 


13, 


.28 


±0. 


.03 


12. 


.51 


±0. 


.03 


11 


.96 


±0. 


.03 


41 


17: 


18: 


:21, 


.22 


-39: 


:18: 


;38, 


.5 


14, 


.86 


±0. 


,04 


13. 


.13 


±0, 


.04 


12 


.01 


±0, 


.05 


42 


17: 


18: 


:25, 


.87 


-39: 


:18: 


:56, 


,9 


13 


.82 


±0, 


,04 


12, 


.71 


±0, 


.03 


12 


.09 


±0, 


.05 


43 


17: 


18: 


:24, 


.70 


-39: 


:18: 


:41, 


,8 


12, 


.71 


±0, 


,04 


12, 


.25 


±0, 


.03 


12 


.03 


±0, 


.05 


44 


17: 


18: 


:25, 


.30 


-39: 


:18: 


:37, 


,8 


13, 


.36 


±0, 


,04 


12, 


.62 


±0, 


.04 


12 


.10 


±0, 


.05 


45 


17: 


18: 


:25, 


.75 


-39: 


:19; 


;25. 


.3 


13, 


.77 


±0. 


.04 


12. 


.81 


±0, 


.04 


12, 


.14 


±0, 


.05 


46 


17: 


18: 


:22, 


.90 


-39: 


:19; 


;40. 


.8 


12 


.91 


±0. 


.04 


12. 


.46 


±0. 


.04 


12 


.18 


±0. 


.05 



Table 2 — Continued 



IRS a(J2000) (5(J2000) JcamlV HcamlV KcamlV 



47 


17: 


18: 


:23, 


.33 


-39: 


:18; 


;27. 


.4 


13 


.09 


±0. 


.04 


12. 


.61 


±0. 


.04 


12. 


.18 


±0. 


.05 


48 


17: 


18: 


:23, 


.35 


-39: 


:18; 


;56. 


.2 


13, 


.32 


±0. 


.04 


12. 


.68 


±0. 


.04 


12. 


.34 


±0. 


.05 


49 


17: 


18: 


:24, 


.48 


-39: 


:18; 


:09, 


.7 


16 


.31 


±0. 


.06 


13, 


.82 


±0, 


.04 


12, 


.51 


±0, 


.06 


50 


17: 


18: 


:26 


.57 


-39: 


19: 


:03, 


,4 


13 


.13 


±0, 


,04 


12, 


,79 


±0, 


,04 


12, 


.55 


±0, 


,06 


51 


17: 


18: 


:24, 


.14 


-39: 


:18: 


:31, 


.0 


13, 


.40 


±0, 


,04 


12, 


,88 


±0, 


,04 


12, 


.62 


±0, 


,06 


52 


17: 


18: 


:32, 


.74 


-39: 


:18: 


:14, 


.0 


12, 


.83 


±0, 


,04 


12, 


,45 


±0, 


.04 


12, 


.33 


±0, 


.05 


53 


17: 


18: 


:23, 


.28 


-39: 


:19: 


;05. 


.9 


15, 


.11 


±0. 


,05 


13, 


.74 


±0, 


.04 


12, 


.63 


±0, 


.06 


54 


17: 


18: 


:21, 


.94 


-39: 


:18; 


;41, 


.4 


14 


.64 


±0. 


,04 


13, 


.37 


±0, 


.04 


12, 


.66 


±0, 


.06 


55 


17: 


18: 


:32, 


.88 


-39: 


19: 


:05, 


.5 


14, 


.56 


±0, 


,04 


13, 


.37 


±0, 


.04 


12, 


.70 


±0, 


.06 


56 


17: 


18: 


:28, 


.32 


-39: 


19: 


:49, 


.8 


14, 


.73 


±0, 


,04 


13, 


.46 


±0, 


.04 


12, 


.71 


±0, 


.06 


57 


17: 


18: 


:23, 


.47 


-39: 


19: 


:52, 


.3 


13, 


.27 


±0, 


,04 


12, 


.93 


±0, 


.04 


12, 


.73 


±0, 


.06 


58 


17: 


18: 


;20, 


.62 


-39: 


19: 


:55. 


.2 


13, 


.59 


±0. 


,04 


13, 


.06 


±0, 


.04 


12, 


.80 


±0, 


.06 


59 


17: 


18: 


:24, 


.74 


-39: 


:19; 


;06. 


.6 


15, 


.00 


±0. 


,05 
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3.2. The IRAS source 

In Figure 6 we show the K band image (CI filter) of the region around IRAS 17149-3916, 
overlapped by the IRAS error ellipse (broken line). Several stars fall inside the ellipse; a more 
accurate coordinate for IRAS 17149-3916 was obtained from the Midcourse Space Experiment (Price 
et al. 2001) - MSX"^ point source catalog. We found one MSX source within the IRAS error ellipse 
(G348.2362-0.9809), with coordinates a(J2000) = 17'"18™25.2", (5(J2000) = -39'^19'^42^ which is 
shown as a contour diagram in Figure 6. None of the stars coincides with the MSX coordinates. 
In fact, from figure 6, we can see that probably several individual sources can be contributing to 
the IRAS emission. 

3.3. Br7 emission from the infrared nebula 

As we already mentioned, there is a strong and compact (with about 1 arcmin radius) radio 
source in the direction of IRAS17149-3916 (PMNJ 1718-3918), with a flux density of 11.9 Jy at 
4.85GHz. In order to study this extended emission, we rc-obscrvcd in the 2005 observing run, the 
field toward IRAS17149-3916 source trough the narrow band K filters, which are centered in the 
Br7 line and in the adjacent continuum (see table 1). 

To calculate the total 617 emission, we first scaled the counts of the continuum image (CI filter) 
to the same level of the Br7 image, using the common bright field stars. In Figure 7, we show the 
Br7 contour diagram constructed from the difference between the Br7 and CI images, overlaying 
the Spitzer^ 3.6/Ltm image of the region around IRAS17149-3916; we also indicated the position of 
four prominent Spitzer sources coincident with the infrared objects IRSl, IRS2, IRS3 and IRS4 
from our survey. Prom that figure, we can also see that the Br7 emission presents a "shell" like 
shape, which seems to be enclosed by 3.6 /im dust emission. From the spectral type of IRSl (05V- 
06V derived from the color-magnitude diagram in figure 5), which is the most luminous source in 
the whole region, we can speculate that the observed Br7 flux density probably is generated by the 



*http://www.ipac.caltech.edu/ipac/msx/msx.html 
^http://www.ipac.caltech.edu/data/SPITZER/SPITZER 



Table 2 — Continued 


IRS 


a(J2000) 


<5(J2000) JcamlV HcamlV 


KcamlV 


93 


17;18;22.1 


-39:19:11.6 14.68 ±0.05 


12.89 ±0.07 


94 


17:18:21.36 


-39:19:54.1 15.60 ±0.06 


13.81 ±0.11 


95 


17:18:30.53 


-39:19:37.6 14.33 ±0.05 


13.82 ±0.11 


96 


17:18:26.59 


-39:17:56.8 13.98 ±0.05 


13.37 ±0.09 



8 



same ionized gas that produces the free-free emission detected at radio wavelenghts. 

Next we constructed the flux cahbrated contour diagrams of the two images (the Br7 and CI 
contour diagrams) and, by measuring the area between contours, we obtained the net Br7 flux 
density using the relation between magnitude and flux density given by Koorneef (1983), finding 
an integrated Br7 flux density of 7.1(=bl.l)xl0~^^ erg cm~^ s~^. 

3.4. Mean visual extinction and emission measure determined from the Br7 and 



To estimate the mean visual extinction Ay in the direction of the nebula, we first used the 
integrated Br7 flux density, uncorrected for reddening, to compute the corresponding number of 
Lyman continuum photons, using the expression derived by Ho (1990): 



We obtained (1.58±0.25) xlO^^ photons s~^. We compared this result with that derived 
from the 6 cm flux density, using the relation from Armand et al. (1996), assuming an electron 
temperature Tg = 7500^: 



radio flux densities 




(1) 
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We found NLy{6 cm) = 5.25 x 10 photons s~^, which corresponds to a mean visual extinction 
Ay = 5.49ibf32 magnitudes. This result agrees with that obtained from the color-magnitude 
diagram for the IRSl source (see fig.4), which appears reddened by about Ay — 6.5 mag. 

We calculate the emission measure E from the detected 5 GHz flux density, using the observed 
size of the Br7 extended emission (about 55 arcsec) and the expression of the expected free-free 
emission from an optically thin plasma at wavelength A, given by: 



= 5.4 X 10-1^5// (A, r)0£;T-i/2g-WAfcT jy 

where E is the emission measure (cm~^), CI the solid angle of the source and gff{X,T) the Gaunt 
factor, which for radio wavelengths can be calculated from: 



5//(A,r) 



V3 



17.7 + In 



r3/2A 



(4) 



Assuming r=7500K and S,y=ll.9 Jy we found £^ = 4.5 x lO^^cm ^, which results in a mean 
electron density rig = 2.6 x lO^cm"^, characteristic of compact HII regions (Churchwell 2002). 



3.5. The IRS4 near infrared source 

In figure 7, we indicated the counterparts of the IRSl, IRS2, IRS3 and IRS4 sources in the 
3.6/im Spitzer image. As we have already mentioned, IRSl is the best candidate to be the main 
ionzing source in the cluster. The other three NIR sources present large excess of emission in the 
color-color diagram, appearing also as very bright IR sources in the 3.6/xm image. 

We attributed the excess emission from IRS2 and IRS3 to dust that probably surrounds these 
young stars. In the case of IRS4, we found a residual in the difference between the point spread 
function photometry (Br7- CI), showing the presence of Br7 line emission, suggesting the presence 
of a highly embedded UCHII region. This can be seen in Figure 8, where we present the K(Br7)- 
K(C1) versus K magnitude diagram, constructed from the point spread function photometry of 
the 41 common bright stars (K <11.5) in the Br7 and CI narrow band images. From that figure we 
can see that all stars but one (IRS4) have compatible K photometry, considering the photometric 
errors. 



4. Conclusions 

NIR observations in the direction of RCW121 reveals the presence of a young cluster of massive 
stars embedded in an HII region coincident with IRAS17149-3916 source. These observations, 
together with published radio data, MSX and Spitzer images were used to determine some of the 
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physical parameters of the region. We found 96 cluster member candidates in an area of about 1.5 
X 2.0 square arcmin, 30% of them showing excess emission in the NIR. IRS 1, the strongest source 
in the cluster with an estimated spectral type of 05V-06V ZAMS based on the color-magnitude 
diagram, is probably the main ionizing source of the HIT region detected at radio wavelengths. 

We did not find a star directly associated with the MSX source coordinates; in fact, probably 
more than one NIR source could be heating the dust, since the strongest of them are concentrated 
near the maximum of the 8.28 fim emission in the MSX image. This result is different from that 
obtained by Roman-Lopes et al. (2003); Roman-Lopes & Abraham (2004a,b), in which a single and 
strong NIR source have coordinates that agree with the MSX source. 

Prom narrow band Br7 and continuum images we measured an extended Br7 flux density of 
(7.1 lb 1.1) xlO"-^^ erg cm~^ s~^. We found that the Br7 emission presents a shell-like structure, 
which seems to be enclosed by 3.6 //m dust emission. Using the integrated Br7 flux density, 
uncorrected for reddening, we computed the corresponding number of Lyman continuum photons 
and compared it with that obtained from the 5 GHz flux density to derive a mean visual extinction 
of Ay = 5.49±f;32 magnitudes. This result agrees with that inferred from the color-magnitude 
diagram for the IRSl source, which appears reddened by about Ay — 6.5 mag. From the observed 
size of the Br7 extended emission and the detected 5 GHz flux density, we calculated emission 
measure E = 4.5 x lO^^cm"^ and electron density n,. = 2.6 x lO^cm"^, characteristic of compact 
HII regions Churchwcll (2002). Moreover, the IRS4 object presents a residual in the difference 
between the point spread function photometry (Br7- CI), showing the presence of Br7 line emission, 
suggesting the existence of a highly embedded UCHII region. 
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Fig. 1. — The J (up), H (center) and K{C1) (down) infrared images (about 5.2' x 5.6'). These 
frames were taken using the near infrared camera coupled to the 0.6 m Boiler & Chivens telescope, 
during the 2003 observing run. North is to the top, east to the left. The "cluster" and "control" 
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Fig. 2. — Color-color diagrams for two comparative regions in our survey. The control diagram 
probably contains only field objects; the nebula diagram has also some field stars but presents many 
young objects with excess emission in the NIR. The locus of the main sequence and the giant branch 
are represented by the solid lines (Koornneef,1983), while the reddening vector for late and early 
type stars are represented by dotted and long dashed lines respectively. A third vector represented 
by the long dotted line indicate the place in the color-color diagram where the stars present excess 
emission in the NIR. 
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Fig. 3. — Comparative histogram of the distribution of the H — K colors in the control (dashed 
line) and in the cluster (solid line) regions. The counts of the control area were normalized to the 
area of the cluster region. 




Fig. 4. — nbK-band contour diagram of the infrared nebula region. The contours are spaced in 
intervals of 10%, 40%, 80% and 100% of the peak value. Here we indicate the position of the 
sources with excess emission (filled triangles), those without (open triangles), the main ionizing 
source candidate (star), and those only detected at H and K bands (filled squares). 



-16- 



8 -r 



10 - 



12 - 



14 - 



16 - 



18 



03v 
04v 
05v 
06v 
07v 
08v 

"09.5V 

BOv 

B0.5V 
B1v 

B2v i 



B1v • 



B5v . 
AOv \ 

A3v i 
A5v ^« 

A8v V 



11 

A 

19 A 



13 

^ 10 



A 7 



\=10 

/ 



8- 



A 



5 

A 



21 T 

27 A A^ 17 15 . 

44 ^ 26 ^■"^ 
Aa34 33 a ^20 
FOv-\ 58^47^483f 31 

68 ''46^43 




A7>"75 65 



A aA 



84 88 



^82 



A 56 
^ 69 



61 

S 55 
57 a60 



54 



A 



'74 

T 

63 



.^'^67 _ 
"^90 64 

T 

77 



▼ 

83 



T 

42 

T 

62 



_L 



2 

T 



■16 

A" 



A 



\ = 20 



22^23. 24 18 

" A. A 



T 

A 

^32^ A- 30 
40 ^ 



12 

A 



A 



36 

A 

50 



▼ 

79 



70 



\ = 30 



-0.5 



0.0 



0.5 



1.0 



1.5 
J-H 



2.0 



2.5 



3.0 



3.5 



Fig. 5. — The J vs J — H color-magnitude diagram of the sources in Tabic 1. The locus of the 
main sequence at 1.6 kpc is shown by the solid line. The intrinsic colors are from Koornnccf (1983), 
while the absolute J magnitudes were calculated from Hanson et al. (1997). The redcning vector 
for a BO ZAMS star {dotted line) is from Fitspatrick (1999). We also indicate the location {bold 
numbers) of Ay = 10, 20 and 30 mag of visual extinction as well as the sources that show excess 
(filled triangles) and those that do not (open triangles) in the color-color diagram. 
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Fig. 7. — Contour Br7 diagram overlayed by the Spitzer 3.6/Lim image of the region around 
IRAS17149-3916. The contours start at 695 counts arcsec"^ and are spaced by the same value. 
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Fig. 8. — The i^(Br7)-i^(Cl) versus K magnitude diagram, constructed from the PSF photometry 
of 41 common bright stars (K <11.5). The anomalous source in this diagram corresponds to the 
IRS4 source (see text). 



